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Summary 

This  paper  considers  analytical  Issues  associated  with  the  notion  of 
the  energy  release  rate  In  quasi-static  elastic  crack  propagation. 


Introduction 

When  a  crack  In  a  deformed  elastic  solid  lengthens,  the  new  free  surface 
created  requires  energy  for  its  formation.  If  the  crack  propagation  process  Is 
slow  enough  to  be  treated  aa  a  quasl-statlc  ,  continuum  mechanical  models  permit 
the  calculation  of  the  rate  at  which  energy  Is  made  available  for  the  generation 
of  new  free  surface  In  terms  of  the  solution  of  the  governing  equilibrium  problem 
corresponding  to  the  deformation  of  the  cracked  solid  under  given  loads. 

The  essential  physical  Ideas  involved  In  this  calculation,  as  well  as  the 
resulting  formula,  have  been  known  for  some  time.  One  may  refer,  for  example,  to 
the  work  of  Atkinson  and  Eshelby  [1],  Budlanaky  and  Rice  [2],  Eshelby  [3,4],  Rice  [5] 
and  Sanders  [6].  To  make  clear  the  most  appropriate  underlying  mathematical  assump¬ 
tions,  however,  seems  to  be  more  troublesome.  Curtin  [7]  has  recently  undertaken 
to  give  a  precise  analysis  In  such  a  way  as  to  be  valid  for  two-dimensional  problems 
In  nonlinear  as  well  as  linear  elasticity.  His  arguments,  however,  rely  on  a 
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The  results  communicated  in  this  paper  were  obtained  in  the  course  of  an  Investiga¬ 
tion  supported  In  part  by  Contract  N00014-C-0196  with  the  Office  of  Naval  Research 
in  Washington,  D.C. 

2 

That  la,  if  inertia  effects  can  ba  neglected. 
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relatlvely  large  number  o£  assumptions,  some  of  which  pertain  to  properties  of 
the  solution  to  the  relevant  boundary  value  problem  which  — ideally,  at  least  — 
should  be  deduced,  rather  than  assumed. 

In  the  present  note,  we  consider  the  simplest  possible  case  of  quasi-static 
elastic  crack  propagation  in  order  to  illustrate  and  clarify  the  purely  analytical 
issues  involved.  We  treat  the  crack  problem  in  anti-plane  shear^  according  to  the 
linearized  theory  of  elasticity.  The  argument  given  here  is  entirely  different 
from,  and  involves  substantially  fewer  assumptions  than,  that  given  in  [7]. 

2.  The  crack  problem 

We  consider  a  homogeneous,  isotropic  elastic  cylinder  in  an  equilibrium 
state  of  anti-plane  shear.  The  outer  boundary  of  a  cross-section  of  the  cylinder 
is  a  piecewise  smooth  simple  closed  curve  C,  while  the  inner  boundary  consists  of 
a  line  segment  of  length  t  (the  crack).  With  cartesian  coordinates  chosen  as  in 
Figure  1,  the  crack  is  described  by  ■  0 ,  .  We  denote  by  R(i)  the 

region  consisting  of  points  on  or  inside  C  which  do  not  lie  on  the  crack. 

According  to  the  linearised  theory  of  elasticity,  the  out-of-plane  displace¬ 
ment  u •  u(x]>x2)  of  particles  in  the  cross-section  of  the  cylinder  is  a  harmonic 
function: 


Au  =  u  +  u  22*^  on  ’ 


(2.1) 


Here  a  comma  followed  by  a  subscript  indicates  partial  differentiation  with  respect 
to  the  corresponding  cartesian  coordinate.  On  the  outer  boundary. 


Anti-plane  shear  refers  to  the  class  of  deformations  of  a  cylinder  in  which  the 
displacement  vector  is  parallel  to  the  generators  of  the  cylinder  and  independent 
of  axial  position.  For  a  discussion  of  anti-plane  shear  (or  "Mode  III"  displace¬ 
ment  fields)  in  the  theory  of  elasticity,  see  [8]. 
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^|n"f(a)  on  C  *  (2,2) 

where  Che  constenc  ^>0  is  Che  shear  modulus  of  Che  elasCic  material,  f 
represencs  Che  given  craccion  on  C,  and  o  is  arc  lengCh ,  measured  posiclve 
counCerclockwlse.  The  derivacive  in  (2.2)  is  in  Che  direccion  of  Che  unic  outward 
normal  vector  n  on  C.  The  faces  of  the  crack  are  to  be  traction-free,  so  that 

u  at  Xj  ■  0±  ,  0  <  x^  <  l  o  (2.3) 


Further, 


u  is  bounded  on  ft(t)  .  (2.4) 

In  the  model  of  quasi-static  crack  propagation  envisaged  here,  the  crack 
lengthens  in  time  by  propagating  to  the  right  only;  i  is  taken  to  be  a  strictly 
increasing  function  of  the  time.  Without  loss  of  generality,  we  may  then  in  fact 
identify  1  with  the  time. 

We  assume  the  existence  of  a  solution  u*u(x;l)  ■ufa^.x^;*.)  to  the 
problem  (2.1) -(2. 4)  which  is  defined  and  three  times  continuously  differentiable 
with  respect  to  (x,l)  on  the  three-dimensional  region  characterized  by 
x€ft(i),  0  <£  *£0  ,  for  some  fixed,  positive  such  that  the  point  (*q.O) 

lies  inside  C.  It  is  further  assumed  that  u  and  all  of  its  derivatives  of  order 
up  to  and  including  three  possess  limits  as  a  point  (x^.O)  on  the  crack  is  ap¬ 
proached  from  above,  provided  Q  <  x^  < t  .  These  limits  are  assumed  to  be  continuous 
functions  of  (x^i)  for  Q<x1<l,  •  Finally,  we  make  analogous  assump¬ 

tions  for  the  corresponding  limits  as  the  crack  is  approached  from  below. 

Except  for  cases  corresponding  to  very  special  choices  of  the  loading  f, 

Vu  is  unbounded  at  both  crack-tips. 
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Ue  shall  first  prove  that  the  energy  In  ft(£)  associated  with  u  , 
defined  by 

E(£)  -  /  ^  |Vu|*dA  ,  0<is*  (2.5) 

R<£)  2  0 

Is  finite  for  each  £  ,  and  satisfies 

2E(£)  -  Jf(o)u(x;£)do  ,  0<£*£Q  (2.6) 

It  follows  immediately  from  (2.6)  that  there  Is  at  most  one  solution  of  the 
prescribed  smoothness  to  the  boundary  value  problem  (2.1) ,  (2.4),  apart  from  an 
arbitrary  additive  constant.  The  boundedness  condition  (2.4)  of  course  plays  a 
central  role  in  the  argument  leading  to  (2.6)  and  the  resulting  uniqueness. 

Our  objective  is  to  study  the  so-called  "energy  release  rate"  £(£)  defined 
by 

£(0 -/f(a)u(x;£)do-E(£)  ,  0<£s»  .  (2.7) 

C 

where  the  dot  indicates  differentiation  with  respect  to  £  .  In  particular 

u(x;£)  «  (x;£)  .  (2.8) 

£(£)  represents  the  excess  of  the  power  supplied  by  the  external  traction  over 
the  rate  of  increase  of  stored  energy,  and  thus  may  be  interpreted  as  the  rate  at 
which  energy  is  made  available  for  the  formation  of  new  free  surface  during  the 
process  of  crack  propagation.  Using  only  the  assumptions  stated  above*,  we  shall 
prove  that 

^ln  particular,  we  make  no  assumptions  about  u  beyond  those  stated  in  the  smoothness 
hypotheses  described  above. 
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e<£)  -  11b  f  (  2  I7u|2,*i  “  fj  u,i^d*  *  (2.9) 

v+0  1^  * 

where  I*  Is  a  circle  of  radius  r  centered  et  the  moving  crack-tip  x^*  t, 

* 0,  s  is  arc  length  on  I*  (increasing  counter-clockwise),  and  n  is  the 
unit  normal  vector  on  I*  which  points  away  from  the  crack-tip;  n^  is  the 

x.  -component  of  n  .  Thus  6(1)  depends  only  on  the  local  behavior  of  Vu 
near  the  moving  crack-tip  (and  not  directly  on  u  ).  We  shall  also  prove  that 
e(l)  i0. 

In  order  to  establish  these  results,  we  need  some  estimates  pertaining  to 
the  behavior  of  u  and  u  near  the  crack-tips.  In  this  connection  it  is  neces¬ 
sary  to  observe  first  that  u  satisfies 

di-0  on  R(£)  ,  (2.10) 

e 

!r-0  on  C  ,  (2.11) 


u  2  •  0  at  x2  •  01  ,  0  <  <  l 


(2.12) 


These  follow  from  (2.1)-(2.3)  by  differentiation  with  respect  to  l  ,  provided 
use  is  made  of  the  fact  that  the  traction  f  is  independent  of  1  .  If  A  were 
bounded  at  both  crack-tips,  and  therefore  on  R(l)  ,  the  uniqueness  result  alluded 
to  above  would  imply  that  O  =  constant  on  S(t)  .  It  will  turn  out,  however,  that 
u  is  in  general  unbounded  near  the  moving  crack-tip. 

We  now  consider  the  local  behavior  of  harmonic  functions  near  a  crack  tip. 
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3.  Crack-tip  estimate  a. 

Let  £  be  the  plane  slit  disc  described  in  polar  coordinates  by 

£■  C(r,0)  |0  <  r  *a,  -  n<  0  <  n}  .  (3.1) 

Further,  let  I*  ,  0<rsa  ,  at  present  be  the  circle  of  radius  a  centered 
at  the  origin,  and  denote  by  ,  0<rSa  ,  and  B^  ,  0  <  r  <  a  ,  the  sets  of 
those  points  in  £  which  lie  respectively  inside  and  outside  r  . 

Let  <p-cp(r,6)  be  a  function  defined  on  £  with  the  following  properties: 

A.  <p€  c\&) ,  and  and  its  derivatives1  cpr  and  cpg  possess  limits  as 
6-*±tt,  0<r*a  ;  these  limits  are  continous  functions  of  r,  0  <  r  . 

B.  (f  is  mean-square  bounded  on  £  ;  i.e.,  there  is  a  constant  m  ^0  such 
that 


2“j"jr  /  <p2ds^m2  ,  0<r^a  , 


(3.2) 


C. 


'prr+~T'pr+  "T  ’86  “  0  on  * 


(3.3) 


D. 


r,±n)  -  0  ,  0  <  r  sa  . 


(3.4) 


We  note  that  <p  has  property  B  if  9  is  bounded  on  £  . 

The  principal  conclusions  concerning  functions  with  properties  A  -  D 
relate  to  the  energy  distribution  associsted  with  cp  and  form  the  basis  of  the 
analysis  to  follow.  The  energy  E(r)  in  £f  associated  with  <p  ,  defined  by 


^Derivatives  with  respect  to  r  or  9  are  indicated  by  subscripts  without 
comas. 
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E(r)  -/  ||Vcji|2dA  ,  0<rsa  , 
r 


(3.5) 


is  finite  and  satisfies 


E(r)  ■  /  7  tr  <fds  20  ,  0  <  r  &a  .  (3.6) 

r  1  r 

r 


Moreover. 


E(r)  -+0  as  r -*-0 


To  prove  this  proposition'1',  we  begin  by  setting 


F(r)  ■  /  #|  V  cp|  2dA  ,  0  <  r  <  a 

«r 


(3.7) 


(3.8) 


A  standard  application  of  the  divergence  theorem,  together  with  (3.3)  and 
(3.4),  leads  to 


F(r)  «g(a)  -g(r)  ,  (3.9) 

where 


g(r)  “  j-  2  cpr'pd8  »  0<ria  (3.10) 

r 


Moreover,  from  (3.9),  (3.8), 


g'  (r)  -  -F'  (r)  -  J  £  |v<p|2ds  20  ,  (3.11) 

r  2 

r 


^The  argument  that  follows  is  essentially  the  same  as  that  used  in  a  more 
general  setting  in  [9] . 
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so  chat  g(r)  is  monotone  nondecreasing.  It  follows  that  - cd  £  g(0+)<  oo 
From  (3.10) , (3.11)  and  Schwarz's  inequality,  one  obtains 


2 

g2(*)  s  4  /  92ds  /  cf,2, 


(3.12) 


By  (3.2)  (property  B) , 


g2(r)  £npm2rg' (r)  ,  0<r£a  .  (3.13) 

Suppose  first  that  g(0+)  <  0  .  Then  for  some  r^  ,  0  <  <  a  ,  one  has  g(r)  <  0 

for  0  <  r  £r^  .  Integrate  (3.13)  from  to  r^,  where  0 < r^ < ,  obtaining 


2  2 
log  —  £  n  pm 

rl 


[8(ri>  '  8<ra)] 


n  pm 


lg(r9) 


(3.14) 


or 


|g(r2)| 


£ 


n  pm2 
logr2/rl 


(3.15) 


Let  r^-*-0  to  8®t  g(r^)  ■  0  ,  which  is  a  contradiction.  Next  suppose  that 

g(0+)>0.  Then  g(r)>0  for  (^rir^  for  some  r^  ,  0  <  r^  <  a  .  Integration 

of  (3.13)  now  gives 


8<Ii> 


(3.16) 


from  which  we  conclude  g(0+)  -  0  ,  again  a  contradiction.  Thus 

g(r)  -*■  0  as  r  +  0  (3.17) 


From  (3.17),  (3.9),  1c  follows  that 

F(0+)  -  g(a) 

is  finite.  But  (see  (3.5)  ,(3.8)  ,(3.9)) , 

E(r)  -  F(0+) -F(r)  -  g(r)  ,  0<rSa  ;  (3.18) 

(3.18),  (3.10)  and  (3.17)  establish  (3.6),  (3.7)  and  the  proof  is  complete. 

4.  Properties  of  u,u  near  the  crack-tips. 

Returning  to  the  original  boundary  value  problem  (2.1)-(2.4),  we  suppose 

first  that  Z  is  a  fixed  but  arbitrary  value  of  Z  with  0 <  Z  <  Zq  .  Let  a  be 

a  constant  such  that  0  <  Z  -  2a  <  Z  +  2a  <  Zq  ,  denote  by  A  the  Z  -  interval 

[Z-  a,  4  + a],  and  note  that  the  (moving)  circle  r*  T  (Z)  of  radius  a  centered  at 

r  r 

the  right  crack-tip  lies  in  the  interior  of  the  boundary  curve  C  and  never  encloses 
the  left  crack-tip,  provided  Z  lies  in  A  .  In  this  section  we  shall  consider  the 
solution  u(x;Z)  and  its  derivative  u(x,i)  only  for  values  of  Z  in  the  interval 
A  . 

With  a  as  above,  let  r^*  1^.(4)  0<rsa  ,  be  the  circle  of  radius  r 
centered  at  x^Z.Xj-O  ,  and  let  &  stand  for  the  set  described  by  (3.1), 
where  r,8  are  now  (moving)  polar  coordinates  centered  at  the  right  crack-tip. 

As  in  the  preceding  section,  denotes  the  set  of  those  points  in  £  which 

lie  inside  ,  0  <  r  <  a  .  Define 

v(r,9;i)  ■  u(Z-4- r  cos  9,  rsin9;i)  ,  0<rSa,-ns9sn  ,  z€  A  ,  (4*1) 


and  note  that 
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v(r . 0; A)  =  gj  (r,6;Z) -u^  +  u  . 


(4.2) 


In  view  of  the  assumptions  made  on  u  in  Section  2,  v  clearly  has  properties 
A-D  for  each  t  in  A  ;  in  particular,  by  (2.4),  there  is  a  constant  m(i)  such 
that 


1 

2nr 


J  v2(r,0;i)ds  *m2(i)  , 

r 

r 


0  <  r  ia  ,  1 6  A 


(4 . 3) 1 


In  view  of  the  proposition  in  Section  3  concerning  functions  with  properties 
A-D,  it  follows  from  (3.7),  (3.5)  with  <p(r,  b)  *  v(r,6;H)  that 


& 


1 7v  j  dA  -*■  0  as  r  -*■  0 


IGA 


(4.4) 


From  (4.2)  and  the  fact  that  u  and  u  are  harmonic  in  &  and  have 
vanishing  first  normal  derivatives  on  the  crack-faces  (see  (2.1),  (2.3),  (2.9), 
(2.11)),  it  follows  that  v  has  properties  A,C,D,  for  each  i  in  A  .  We  show 
next  that  v  has  property  B  (see  (3.2))  for  each  such  l  .  To  this  end,  let 
l  and  l'  belong  to  A  ,  with  %'  +  l  .  Define 


v(r>6;i»i,)“ 


v(r,0;JO-v(r.8;i/ ) 
l-  i' 


0<ria,  -tis0sn,  l,l' €  A,  l'  i  l  .  (4.5) 


Clearly  V  has  properties  A-D,  and  hence,  choosing  cp=V  in  (3.6)  and  making 
use  of  (3.7),  one  has 


JwdsiO,  0<rsa  ,  &,«.'€ A  ,  l  + 1'  .  (4.6) 

r  r 

r 


The  idea  now  is  to  let  l'  ■*  l  with  r  fixed  in  (4.6).  To  carry  out  this  limit 


^In  fact,  j v |  *  m  on  £  . 
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process,  we  first  observe  from  (4.5)  that 

V(r,0;i,i  )  -*-v  r,0;£)  as  £■+•£,  for  each  fixed  (r,9)  .  (4.7) 

Moreover,  from  (4.5),  the  mean  value  theorem,  and  the  smoothness  properties  that 
v  inherits  from  u  through  (4.1),  it  follows  that 

|v(r,9;£,£')!  SM^r)  ,  0<rSa  ,  -nsfisn  ,  £',£€A,  l' + 1  .  (4  g) 

and 


fr(r,0;£,£')|  ^M2(r)  ,  0<r*a,  -ns9sn  ,  £',£6  A,  £'  +  £ 


(4.9) 


where 


ML(r) 


max  |  v(r,0;  X)  |  ,  M„(r)  »  max  |v  (r,0;X)|  ,  0<r 

-tsesu  z  -ts0sn 

X€ A  16  A 


sa 


(4.10) 


It  follows  that 

|w  |  on  for  each  r,  0<rSa,  and  for  all  £*,£€  A,  l  f  £  .  (4.11) 

Thus  the  polntwise  convergence  in  (4.7),  although  not  necessarily  uniform,  is 
dominated.  It  follows  from  a  standard  theorem^-  that  the  limit  as  l'  -*■  l  may  be 
taken  inside  the  integral  in  (4.6),  yielding 


/vvdsiO,  0<rsa  ,  £6  A  . 


(4.12) 


Thus 


X£ee  [10],  p.  405 


-12- 


0  <  r  ia 


i€A  , 


(4.13) 


so  that  the  integral  in  (4.13)  is  a  nonnegative,  monotone  nondecreasing  function 
of  r  .  It  follows  that  this  integral  is  a  bounded  function  of  r,  0  <  r  &a  ; 
i.e.  there  is  a  constant  n(£)  *0  such  that 


TSF  /i2d«‘n2U)  , 

r 


0<rSa  ,  A  . 


(4.14) 


Thus  v  has  property  B  for  each  l  in  A.  The  main  proposition  of  Section  3 
may  now  be  invoked  with  <f-v  ;  from  (3.6),  (3.7)  we  find  that 


/  |Vv|2dA  +  0  as  r  +  0  ,  i€A  . 

We  note  that  (4.3),  (4.14)  imply 


/  v2ds  -  0(r) 

r 

r 


/  $  ds  ■  0(r)  as  v-^0  ,  i.€A  . 

r 

r 


(4.15) 


(4.16) 


Finally,  we  observe  that,  since  the  choice  of  the  center  i  of  the 
Interval  A  was  arbitrary,  each  of  the  four  limiting  results  (4.4),  (4.15), 
(4.16)  holds  for  every  l  in  (0, 1^) . 

It  may  be  remarked  that  a  more  refined  argument,  using  more  detailed 
properties  of  v  ,  can  be  used  to  prove  that  $  is  actually  continuous  at 
the  moving  crack-tip11,  and  not  merely  mean-square  bounded  as  asserted  in  (4.14). 

^Continuity  of  the  analogous  physical  quantity  in  plane  strain  at  the  moving 
crack-tip  is  assumed  in  [7]. 


Mfe***.,  • 


f 


Since  (4.14)  is  sufficient  for  our  purposes,  we  omit  the  proof  of  the  continuity 
of  v  . 


As  remarked  earlier,  Vu  is  in  general  not  bounded  —  or  even  mean-square 
bounded  —  near  a  crack-tip.  According  to  (4.2)  and  (4.14),  therefore,  u  would 
in  general  not  be  bounded  near  the  moving  crack-tip. 

Analysis  parallel  to  that  given  above  can  be  carried  out  near  the  left 
(fixed)  crack-tip.  One  sets 

w(r,  0;i)  •  u(r  cos  0,r  sin  b;i)  ,  (4.17) 

where  r,0.  are  now  polar  coordinates  centered  at  the  origin.  In  contrast  to 
(4.2) ,  we  have 

w*u  (4.18) 

If  r.  is  a  circle  of  radius  r  centered  at  the  origin  while  is  the 

interior  of  with  the  crack  deleted,  one  can  again  show  that 


/  w  ds  "0(r) 


/  wds»0(r)  as  r  +  0  ,  0  <  l  <  tn  , 

r  ° 

r 


as  well  as 

2  2 

/  |  Vw|  dA  -*  0  ,  /  1 7w|  *dA-*-  0  as  r  -*0  ,  0  <  l  <  Zn 

A  A  ^ 


He  omit  the  details. 


(4.19) 


(4.20) 


^ Again,  w  can  be  shown  to  be  continuous,  and  not  merely  mean-square  bounded, 
at  the  left  crack-tip. 


IC  is  clear  from  Che  local  considerations  of  Che  preceding  section  that 

the  energy  E(£)  of  (2.S)  associated  with  Che  boundary  value  problem  (2.1)-(2.4) 

is  finite.  In  order  to  study  E (!)  ,  we  let  be  circles  of  radius 

r  ,  0<rSa  ,  centered  at  the  left  and  right  crack-tips,  respectively1.  Let 

Rr(l)  be  the  set  of  those  points  which  lie  in  R(i)  but  outside  , 

0<rsa  .  An  application  of  the  divergence  theorem  to  the  integral  of 
2 

( p/2)  |  Vu|  over  Ml)  ,  accounting  for  the  presence  of  the  crack  and  making 
use  of  (2.1),  (2.3),  followed  by  the  limit  r-*>0  ,  leads  easily  to  the  conclusion 
that 

E(t)  -  \  I  f(o)u(x;A)da  ,  0<UL  ,  (5.1) 

z  C 

as  claimed  in  (2.6). 

Since  the  given  traction  f  on  C  is  independent  of  l  ,  it  follows 
immediately  from  (5.1)  that 

2E(Z)  ■  /  f(c)u(x,i)do  .  (5.2) 

C 

An  alternative  expression  for  E(&)  may  be  obtained  from  (5.2)  by  proceeding  as 
follows.  From  (5.2),  (2.2), 


Tiers  a  is  assumed  small  enough  to  keep  and 

assure  that  each  circle  encloses  only  one  crack  tip. 


inside  C  and 


to 
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2E(1)  */l*|^udo-/  u|^ud8 


3n 


t  3u  .  , 

-I 


di) 

r 


-  I  V>  u  ds  , 

r(2)  3n 


(5.3) 


where  S  ■C+I'^  +  ,  and  n  is  the  unit  normal  on  S  ,  taken  outward 

r  r  r  ***  r 

with  respect  to  <1^(1).  The  arc  lengths  and  integrals  over  ,  1^^  are 

taken  counterclockwise.  The  divergence  theorem,  together  with  the  facts  that 
u  and  u  are  harmonic  and  satisfy  the  free  surface  conditions  (2.3)  and  (2.12) 
on  the  crack  faces,  permits  one  to  write 

/p~uds-/|*-~uds  .  (5.4) 

Sr  Sr 

Making  use  of  (2.11)  in  (5.4),  we  find  that 

/  y.|^uds-J  »i|£uds+J  n|£uds,  (5.5) 

Sr  3n  jU)  3,1  j<2)  3n 

r  r 


so  that  (5.3)  may  be  written  in  the  form 

2E(i)  -I^rji)  +I2(r,i) 


(5.6)1 


where 


*The  representation  (5.6),  (5.7)  for  E(t)  is  the  analog  in  anti-plane 

shear  of  Sanders'  I-integral  in  plane  strain  or  plane  stress  (see  Eq.(4)  of  [6]). 


[»j< 


Equations  (5.6),  (5.7)  show  that  E(£)  depends  only  on  the  behavior  of  u 
and  u  near  the  crack- tips . 

Consider  first  the  integral  :  from  (4.1)  and  (4.2), 


V2J  +  H2  , 


where 


J- J(r 


'‘’'{.(a)  2  [lS",l  ds  > 


H,  «H.(r;i)  ■  /  p<vv-vv)ds. 
22  ^2)  r  r 
r 


(5.8) 


(5.9) 


(5.10 


A  direct  calculation  which  makes  use  of  the  facts  that  v,  v  are  harmonic  and 
vg,  Vg  vanish  at  9«+tj  reveals  that  dl^/fc^O,  so  that  the  integral  in 
(5.10)  is  Independent  of  rsH^r;!)  -l^d).  If  one  now  integrates  (5.10)  with 
respect  to  r  over  the  interval  0  <  r  <  r^  ,  r^  <  a  ,  one  finds. 


r^d)-/  p(vrv-vrv)dA 


The  Schwarz  inequality  now  gives 


r?H»(i)  *2^  {/  v2dA  /  v*dA  +  /  v2dA  /  v2dA  ] 

1  ^  A  *  ^  ^  *  A 


(5.11 


(5.12 


*Tha  finiteness  of  the  Integrals  on  the  right  in  (5.12)  assures  the  existence  of 
the  integral  in  (5.11). 
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Froo  (4.14)  and  (4.3)  one  finds  Chat 


f  v^dASnn^(l)t? 


so  chac  (5.12)  yields 

H2(i)  fm2(l)/  | VO 1 2dA  +  n2 ( A) /  |?v|2dA 

L  &  4 

rl  rl 

Letting  r^-*0  in  (5.14)  and  using  (4.15),  (4.4),  we  find 

h2U)  -0  , 


(5.13) 


(5.14) 


(5.15) 


so  that,  from  (5.8)  , 


lim  I  (r;i)»lim  J(r;i)  (5.16) 

r-*>0  2  r-*-0 


As  to  the  integral  I2(r,l)  ,  one  shows  by  a  similar  argument  based  on 
(4.18),  (4.19)  and  (4.20)  that 


11m  I. (r;i)  »  0 
r  -*■  0  1 


(5.17) 


Letting  r  +  0  in  (5.6)  and  using  (5.16),  (5.17),  and  (5.8),  we  obtain 


E(i)  -  lim  J(r;A) 
r  -*-0 


(5.18) 


where  J  is  given  by  (5.9).  Note  that  J  does  not  Involve  &  and  depends  only 
on  the  behavior  of  u  near  the  moving  crack-tip. 

A  more  convenient  representation  for  J  can  be  found  as  follows.  Since 


(4.1), 


•18— 


u  ,  ■  v  cos  0  -  — 
,1  r  r 


sin  6  , 


(5.19) 


one  nay  write  (5.9)  as 


“  2 


/ 

^ 2 ) 


[  v  (v  cos  9 - 

rv  r  r 


/gain  0)-  v  ~  (vf  cos  6  -  —  VgSin  9)  ]  ds  (5.20) 


The  fact  that  v  is  harmonic  can  be  used  to  eliminate  v  from  (5.20);  integra¬ 
tions  by  parts  and  the  boundary  conditions  ■  0  at  0>±tt  then  give 


J-J 

r<2) 


Oil  V  . 

“,i  jT>  d* 


(5.21) 


Here  n  is  the  unit  normal  on  P  '  pointing  away  from  the  crack-tip,  and 
n^  ■  cos  0  is  its  -  component .  The  integral  in  (5.21)  is  the  usual  J-integral 
of  fracture  mechanics  discovered  by  Eshelby  [4]  and  Independently  by  Rice  [11]. 

As  these  authors  have  shown,  and  as  is  easily  verified,  J  is  independent  of  r. 
In  fact,  the  integral 


J- J(i)  -  J  p(  j  IVul^-u  ^  )  ds  ,  (5.22) 

is  Independent  of  the  path  T  ,  as  long  as  T  is  a  curve  in  frf.1)  which  begins 
and  ends  on  the  faces  of  the  crack  and  surrounds  the  right  crack-tip  but  not  the 
left  one.  Thus  (5.18)  may  be  replaced  by 

fi(t)-JU)  ,  (5.23) 

with  J(t)  given  by  (5.22).  In  view  of  (5.2),  (5.23)  the  energy  release  rate  of 
(2.7) is  given  by 
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e(D-fcU)-JU)  (5.24) 

Thus  Che  power  delivered  by  the  external  load  is  divided  equally  between  the 
casks  of  increasing  the  stored  energy  and  making  energy  available  for  crack 
propagation. 

Equations  (5.24),  (5.22)  show  that  £(i)  (and  E(i))  depend  only  on  the 
local  behavior  of  7u  near  the  moving  crack-tip. 

In  order  to  show  that  £(£)  20  ,  it  is  convenient  to  introduce  a  function 
&(l, l')  as  follows:  for  0<i/si<iQ  ,  set 

<?<M')«J  $  |7u(x,i')|2dA- Jf(c)u(x,i')  do  ,  (5.25) 

8(i)  C 


and  let 


PU)  -  $U,i)  ;  (5.26) 

P(i)  is  called  the  potential  energy1  in  R(i)  associated  with  u.  Moreover, 
(5.1),  (2.5)  give 


P(i)--EU)  , 


(5.27) 


so  that 


PU)--E(i)  . 


(5.28) 


Let  if)  ,  if  >  be  circles  of  radius  r  centered  at  (0,0)  and 
(£,0)  as  before,  and  let  if)  be  a  circle  of  radius  r  centered  at  (*,'  ,0)  , 

It  follows  from  properties  of  uOc,£)  already  established  that  ^(1,1*),  P(l) 
are  finite 
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whsre  0<l'<l  .  Let  8>r(i,t')  be  the  Mt  of  those  points  in  8(1)  which 

lis  out  aid*  r*1*  +  r*2)  +  I*3)  .  Sat 
r  r  r 

*>(£,£')-/  |  l?(x,l')|2dA-/  f(o)u(x,l')do  ,  0<l'  <!<!.,.  (5.29) 

r  «ra.t)J  C 

The  divergence  theoren  and  (2.1)  -  (2.3)  sake  it  possible  to  show  that 

*  (l,l')-*  (1,1 )-/  ,  #  |7u(x,l)-7u(x,l')|2dA 

r  r  Rr(l,l  )  2 

+  Qr(t,l')  *Qr(i,i')  ,  (5.30) 

where 

Qjl.l')-  I  fu(*,t') -u(j.,l)J  Qc,£)  ds  .  (5.31) 

’  rU>+r(2)+r<3) 

r  r  r 

By  using  the  crack-tip  estlaates  of  Section  4,  one  can  show  that  Qr(l,£#)-*-0 
as  r-» 0  .  Letting  r-*-0  in  (5.30)  and  noting  ^  (i,l#)  4^(1, 1#)  as  r-*-0, 

we  obtain 

9  (!,*')  *9  (1,1)  -P(l)  ,  0<l'<£.  (5.32) 


But 

/  |  |Vufe,l')|2dA-/  jl^uQc.l' )  |2dA  , 

8(1)  8(1 ') 

so  fron  (5.25),  (5.26),  (5.33), 

9  (l,l')-f  (t',i}-P(l')  ,  0<l'<l  . 


(5.33) 


(5.34) 


f 


Combining  (5.34) ,(5.32)  gives 

P(i')2P(i)  0  <  i#<  l  ,  (5.35) 

so  Chat  P( £)  is  nonincreasing  with  l  .  Thus  from  (5.28),  E(i)  20  ,  and 
from  (5.24),  &(l)  20  . 

The  inequality  (5.35)  may  be  viewed  as  an  instance  of  the  minimum  principle 
associated  with  the  boundary  value  problem  (2.1)-(2.4)  (the  principle  of  minimum 
potential  energy) . 
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